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Abstract. We consider the dispersion correction to elastic parity violating electron-proton scatter- 
ing due to yZ exchange. In a recent pubUcation, this correction was reported to be substantially 
larger than the previous estimates. In this paper, we study the dispersion correction in greater detail. 
We confirm the size of the disperion correction to be ^6% for the QWEAK experiment designed 
^ to measure the proton weak charge. We enumerate parameters that have to be constrained to better 

than relative 30% in order to keep the theoretical uncertainty for QWEAK under control. 
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I Precision tests of the Standard Model at low energies provide an important frame- 

^ work for New Physics searches and for setting stringent constraints on the parameters 

4^ of possible extensions of the Stadard Model (SM). Such tests involve high precision 

measurements of parameters that are suppressed or precisely vanish in SM. An impor- 
^ tant example of such parameter is the weak charge of the proton, Q^^ = 1 — 4sin^ Ow- 

^ With the value of the weak mixing angle at low momentum transfers sm^Ow{0) = 

O 0.23807 ±0.00017 [1], the SM predicts the proton weak charge of order ^ 0.05. A 

^ precise 4% (combined 2% experimental and 2% theoretical uncertainties) determination 

of the weak charge of the proton is the aim of the QWEAK experiment at Jefferson 
2 Lab [2]. In this experiment, parity violating asymmetry in polarized electron scattering 

O ^pv = — Ol) / {(^R + Ol) will be measured, Or{ol) standing for the differential cross 

^ section in elastic electron-proton scattering with the incidental electron beam spin paral- 

> lei (antiparallel) to its direction, respectively. To leading order in the momentum transfer 

^ t and in Fermi constant, the parity violating asymmetry arises from the interference of 

H parity conserving and parity violating amplitudes at tree level. Including radiative cor- 

rections of order a, one has 

= ^^^!2^[1 +Re5^c + Re5^(v)] + ^(r^) , (1) 

Above, we denote all the radiative corrections except the yZ direct and crossed boxes 
by 5rc, and we explicitly indicate that the dispersion correction dyz is function of energy 
V. The corrections 5rc were considered in various works, to mention the most important 
references [1, 3], and the combined theoretical uncertainty associated with these was 
shown not to exceed 2.2%. The dispersion correction dyz is the only one that can obtain 
a sizeable contribution from hadronic structure. In [4], the dispersion correction was 



represented as the sum 5^ = b^z^ + ^yZv> and two forward sum rules were established. 
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The respective imaginary parts can be expressed as integrals over the PV DIS structure 
functions Fi{x, Q^), 
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with = {M + mnY the pion production threshold, and Q^^ax = (-^ ^ K-^ ^ ) _ 
sum rules of Eqs. (2,3) are new results. In the past, only their values at v = were 
calculated in the framework of atomic PV. Due to the explicit factor of V, the correction 
5yza is exactly zero for zero energy, in accordance with [3], and has been neglected 
in the literature even for non-zero energies. The dispersion relation calculation of [4] 
lead to dyz^ ~ 6% that should be compared with zero in the analysis of the QWEAK 
experiment. The correction dyzy obtains its value due to hard kinematics inside the loop 
and is largely energy-independent [3, 5]. In the rest of this article, we concentrate on 
5yza only. The sum rule of Eq. (2) itself is model-independent; however, in absence of 
any detailed PVDIS data, the input in this sum rule will depend on a model. We proceed 
by modeling the electromagnetic data first. 



(4) 



MODELING REAL AND VIRTUAL PHOTOABSORPTION DATA 

We will use the following three models to model the electromagnetic DIS structure 
functions Fi: 

• Model I: The model used in [4] utilized the resonance parameters obtained in [6] 
and the non-resonant Regge contribution from [7] that was fitted to the real photon 
data at high energies. Ref. [7] also provides the Q^.^ependence of the high-energy 
part, so no additional modelling was necessary here. For the estimates of [4], a 
simple dipole model with the dipole mass A 1 GeV for all the transition resonance 
form factors was employed. 

• Model II: Another form of resonance and background contributions and transition 
form factors was used in [15] to fit virtual photon data. 

• Model III: As Fig. 1 clearly demonstrates, the background of [15] cannot be used 
too far beyond the resonance region. Therefore, we opt for Model III that uses the 
background from Model I, and adopt the resonance contribution from Model II. 
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FIGURE 1. World data on total photoabsorption [8, 9, 10, 11, 12, 13] (see [14] for the complete list) 
compared to the three models described in the text. 



In Figs. 1 and 2, we display the comparison of the three models with the data for the 
differential cross section for inclusive electroproduction in the resonance region. The 
Model I [4] is shown by solid red lines, IVIodel II by the solid blue line, and IVIodel 
III is represented by a long-dash red line. From Fig. 2 it becomes clear that IVIodel I 
largely underestimates the JLab data in the whole resonance region, apart from the real 
photon point Fig. 1 . IVIodel III also provides a very good description of the DIS data for 
Fiix, Q^) at low and moderate values of [7]. 

ISOSPIN STRUCTURE 

Once the parameters are fixed, we can employ simple isospin considerations for each 
contribution. Note that because breaking down the experimental data into the sum 
of resonances and background is a model-dependent procedure, also the isospin ro- 
tation is only defined within a given model. For the N N*{I = 1/2) transition, 
the isospin decomposition resembles that for the elastic form factors, {N*\JI^(-- y\p) = 

(1 — 4s^6w){N*\Jem\p) — {N*\Jem\n) It is then straightforward to relate the contribution 
of a resonance R with isospin 1/2 to the interference yZ "cross section", to its contribu- 
tion to the electromagnetic cross section by introducing the isospin scaling factors 

yZ,p aP ah* iaP An* 

pR _ ^T,R .2n ^ ^R.M2^RM2^^R,3/2^R,3/2 

^z/r=-znp = (1-45%) -p — r^—p — ^ — (5) 

Above, ^^^"/2(3/2) transition helicity amplitudes for exciting the resonance 

R on the proton (neutron), respectively, whereas Ot^r denotes the contribution of the 
resonance R to the transverse virtual photon cross section. The scaling factors 





FIGURE 2. Differential cross section data in the resonance region from JLab E002 [16] are shown in 
comparison with the three models. 



depend both on the relative size and relative phase of the transition amplitude on the 
proton and on the neutron. To obtain both, we use the results of the constituent quark 
model of Ref. [18] for the transition helicity amplitudes. For isovector transitions A'^ — )■ 
A(A*), the scaling factor in the SM is i^^/y = 2(1 — 2 sin^ 6w)- Vector Meson Dominance 
Model (VDM) capitalizes on the fact that the photon (Z-boson) has the same quantum 
numbers as vector mesons and can be represented as a superposition of a few vector 
mesons, I7) = Y.v=p.a),(j>,...\^)- the naive VDM, these three channels cannot mix 
among each other, and one obtains a prediction for the ratios of the cross sections of 
p, (0, (^) production, (J^^^^PP : o^p^'^p : o^P^fP = l : {q'=^/q'=^)^ : {q'/q'=^Y = 
1:1/9:2/9, where we defined q'=^ = 1 /(3v^), q'=^ = 1 /V2 and q' = -1/3. The 
above predictions were confronted to the experimental data at high energies and for 
that ranged from zero to several GeV^ [19] and showed a very good agreement for the 
(o/p ratio, while for the 0/p ratio the agreement was not as good. Accomodating these 
considerations here, we obtain the following ratio of the high energy non-resonant (NR) 
contributions to 7*p — t- Zp and y*p ^ y*p cross sections: 



^^ly ^ (^/=i)2 + (^/=0)2 + (^.^2 -2(l-4sm e^) (6) 



where we make use of the SM definition gy'^ = — 4sin^ 0iy/(3-\/2), q^ ^ = (2 — 
4 sin ew)/V2 and ^-^ = - 1 + 4 sin^ dw /3. We summarize these results inTable 1 . 



TABLE 1. Isospin scaling factors ^^/y for resonances and non-resonant (NR) background. 



P33(1232) 


5ii(1535) 


£'13(1520) 


5ii(1665) 


i^i5(1680) 


^11 (1440) 


F37(1950) 


NR 


1.075 


0.885 


0.938 


0.473 


0.35 


0.745 


1.075 


1.075 




FIGURE 3. Results for the dispersion correction as calculated in the Models I, II and III. 



RESULTS FOR Redyz^ 

Fig. 3 shows the energy dependence of Re Syz^, calculated in three different models. The 
individual contributions at QWEAK energy are shown in Table 2. We next discuss the 
size and origin of the uncertainty in calculating Re Syz^ ■ It comes about from i) modeling 
the electromagnetic data, and ii) from the respective isospin structure. To estimate the 
former, we average over the three models discussed above and displayed in Fig. 3, 
leading to ReSyz^ = (5.85 ±0.45)%, with most of the error bar due to the Model 1. 
Because the latter underestimates the experimental data throughout the resonance region 
for non-zero Q^, it is clear that this estimate represents an upper bound for the model 
dependence. In an upcoming work, we plan to directly fit the parameters of Model III 
to the data that will allow to significantly reduce the model-dependent error, given the 
quality of the experimental data shown in Figs. 1, 2. The isospin structure represents 
a more significant source of uncertainty. We will next indicate, to what precision the 
isospin scaling factors of Table 1 should be constrained to keep the uncertainty of Re dyZA 
under control. From Table 2, it is seen that contributions of Di3(1520), Fi5(1680), 



TABLE 2. Individual resonances and background (NR) contributions to the dispersion correction 
ReSyZA for QWEAK energy Eiat = 1.165 GeV. Results of Model III are quoted. 



^33(1232) 


5ii(1535) 


Di3(1520) 


Fi5(1680) 


5ii(1665) 


Pi 1 (1440) 


^37(1950) 


NR 


1.27 % 


0.44 % 


0.21 % 


0.06 % 


0.05 % 


0.10% 


0.50 % 


3.61 % 



5*1 1 (1665) and Pi i( 1440) resonances have very little impact on the size of the dispersion 
correction. We can therefore afford to assign a conservative uncertainty of 50% to each. 
^33(1232) has been studied in many phenomenological models, and we believe that it is 

reaUstic to constrain the scaling factor ^^/^^^^^^ to better than 20% that would translate 

into 0.25% for the dispersion correction. For 511(1535), the isospin structure might be 

more complex, since it couples strongly to kkN states along with kN, but the absolute 
size of the correction quoted in Table 2 allows for a larger relative uncertainty here. 

Constraining i^^'y to relative 30% would lead to 0.13% for 5yz- Finally, we turn to 
the uncertainties associated with the higher energy contributions. It has to be mentioned 
that the contribution associated in [15] with a resonance around W = 1950MeV cannot 
be ultimately identified with the ^37(1950) state listed in PDG [14]. This contribution 
may represent a missing strength of the background at W > 1900MeV that can be 
seen in the two lower panels of Fig. 2, rather than a real resonance. Constraining this 
contribution to relative 50% means 0.25% for 5yz- For the background, we assign the 
100% uncertainty to the strange quarks contribution to the scaling factor in Eq. (6), 
1^^^ = 1.075(1 ±0.34). Putting the above discussion together, we write 

Re%(v= 1.165 GeV) = (5.85±0.45^orf± l-23;.«)% • O) 

We conclude that the main uncertainty to the dispersion correction comes from the 
isospin decomposition of the electromagnetic data, most notably from high energy back- 
ground. We showed how the total uncertainty can be represented as an uncertainty in the 

four isospin scaling parameters, ^^/y^ ''z/y ^z/y^^^^^ ^^'^ '§2/7'^^^^*^^ " estimates 
of uncertainties in these factors that we presented are only approximate, and will be 
studied in greater detail in an upcoming work in order to ensure the interpretability of 
the QWEAK experiment. 
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